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INTRODUCTION 

For  designing  extractive  or  azeotropic  distillation  systems 
a  knowledge  of  ternary  vapor-liquid  equilibrium  data  is  required. 
There  is  a  large  number  of  binary  vapor-liquid  equilibrium  data 
in  the  literature,  but  few  ternary  systems.  The  experimental 
determination  of  vapor-liquid  equilibrium  data  for  ternary  systems 
requires  long  and  detailed  experimental  work,  and  is  always  com- 
plicated by  inaccuracies  in  the  methods  of  analysis. 

Therefore  a  simple  and  systematic  method  of  calculating  ac- 
tivity coefficients  in  ternary  systems  from  binary  equilibrium 
data  available  for  the  three  binary  combinations  of  components 
would  be  quite  useful.  Such  a  method  would  make  possible  the  re- 
liable prediction  of  vapor-liquid  equilibrium  data  for  any  ternary 
system  with  the  aid  of  binary  vapor-liquid  equilibrium  data  alone, 
or  at  worst,  only  a  relatively  small  number  of  experimental  mea- 
surements on  the  ternary  system. 

Several  authors  have  attempted  to  correlate  vapor-liquid 
equilibrium  data  of  non-ideal  solutions  on  a  semi -empirical  basis, 
notable  among  them  being  Margules  (10),  van  Laar  (24),  Scatchard 
and  Hamer  (15),  Carlson  and  Colburn  (3),  White  (26),  Benedict 
et.  al.  (1),  and  Wohl  (27).  Wohl  combined  all  these  authors' 
attempts  into  a  single,  general  equation. 

This  work  was  concerned  mainly  with  Wohl's  Margules  three- 
suffix  equation  and  whether  ternary  activity  coefficients  could  be 
predicted  from  binary  data  alone,  or  whether  a  ternary  constant  was 


necessary.  Non-aqueous  systems,  taken  from  the  literature,  were 
used  in  this  investigation. 

The  Margules  binary  constants  were  determined  by  several 
methods  and  every  combination  was  used  for  the  ternary  calculations, 
The  activity  coefficients  of  each  component,  which  are  calculated 
by  Wbhl's  Margules  equation  with  various  binary  constants,  were 
compared  with  binary  and  ternary  experimental  data  in  order  to 
determine  which  method  for  the  calculation  of  binary  constants  is 
the  best. 

Ternary  constants  were  thought  to  be  necessary  whenever  the 

average  of  the  absolute  value  of  the  deviations  of  the  calculated 

activity  coefficient  from  the  experimental  activity  coefficient 

was  more  than  or  equal  to  ten  percent.  (  e.g.,  4*-  t&   i  -"/ei)>0>l   ) 

isi  ntfei 

MATHEMATICAL  REPRESENTATION  OF  ACTIVITY  COEFFICIENT 

The  activity  coefficient  "o  ±   may  be  thought  of  as  the  ratio 
between  the  apparent  vapor- pressure  of  the  ith  component  in  a  li- 
quid solution  to  its  Raoult's  law  value.   In  an  ideal  solution,  the 
activity  coefficient  of  every  component  equals  one,  and  in  any 
type  of  solution  its  value  for  a  particular  component  approaches 
one  as  that  component's  mole  fraction  approaches  one.   In  non- 
ideal  solutions  activity  coefficients  are  strong  function  of 
concentration,  but  their  variation  must  satisfy  the  Gibbs-Duhem 
equation  for  thermodynamic  consistency. 

In  the  preceding  section  it  was  mentioned  that  Wohl  generated 


several  authors'  (1,3,  10,  15,  24,  26,  27  )  equations  from  a 
single  general  expression  relating  the  activity  coefficients  as 
functions  of  concentration  in  binary,  ternary,  and  quaternary 
systems.  Using  the  excess  free  energy, A ge,  as  a  function  of 
composition,  he  obtained  equations  suitable  for  expressing  the 
activity  coefficients  of  the  components  as  a  function  of  the  com- 
position of  the  mixture, 

Vfohl  expressed Aqe  in  terms  of  the  following  expansion; 


S03 

where 


2.3d3RT£qixi  =4ZiZjaij   ^Vj^ijk  +^,ZiZJZkZlaiJkl    (2-1) 


JJ  AJK  i.jk£ 


&G,     is  molar  excess  free  energy; 

3i»  3j  are  the  effective  volumes  of  the  components  i,  j; 

aij»  aijk»  ^jkl...  are  empirical  constants  measuring  the  in- 
teractions in  various  groups  of  molecules  ij,  ijk,  ijkl 

The  effective  volumetric  fraction  of  any  component  i  is  defined 
by  the  relation, 

zi  =       qixi  (0  9s 

23ixi 
Equation  (2-1)  is  written  as  an  equation  of  fourth  order. 

From  this  equation  Wohl's  activity  coefficient  equation  of  fourth 

order  for  component  k  can  be  obtained.   It  is  usually  found  that 

the  higher  the  order  of  the  equation  used,  the  better  it  will 

represent  the  behavior  of  the  given  system,  but  the  more  complicated 

it  will  become.  Practically  speaking,  the  equations  of  order 

higher  than  four  are  rarely  used,  and  the  equations  of  third  order 

were  used  in  this  work. 


(A)  Binary  Systems.  An  equation  of  the  third  order  of  the 
form  of  equation  (2-1)  can  be  written  for  a  binary  system, 

2^&R't  =   <   xl   +  ^-  x2    )Z2Zl  [ziql(2ai2   +  3an2) 

+  Z2q1(2aI2  +  3a122)J*  <2-3) 

After  introducing  the  new  constants, 

A12  a  qi (2*1.2  +  3a122) 

(2-4) 

*21  =  ^2^*12  +  3a112) 
equation  (2-3)  can  be  written  as: 

^g^  -      (   X!   ♦  32x2    )Z2ZL(A2lZifjl  +  Ai2Z2    ).  (2-5) 

Wohl's  activity  coefficient  equation  of  the  third  order  for  a 
binary  system  can  be  obtained  from  equation  (2-5); 

logflV  =  Z22  Ca12  *  2Z1(A21qi/q2  -  Ai2)J, 

log^2  =  Z^  |^A21  +  2Z2(A12q2/qL  -  A21)] 

On  the  assumption  that  q2/qisslt  equations  (2-6)  take  the  form 
of  the  three-suffix  Margules  equations.  These  equations  contain 
two  constants  A^2,  A2^  which  have  to  be  determined  from  the  ex- 
perimental measurements: 

logl^  ■  x22  £a12  +  2xL(A21  -  A^))  , 

log\2  =  X],2  j^A2l  +  2x2(Al2  -  A21)|  . 

Redlich-Kister  equations  (13)  are  rearranged  forms  of  the 
three- suffix  Margules  equations: 

logT^  ■  x22  [b12  +  c12(x2  +  Sx^J  , 

logJ2  =  ^2  [b12   +  c12(xL  -  3x2)J  . 
The  relations  between  Redlich-Kister  and  Margules  constants  are 

b12  *  £  (A12  +  A21),  (2-9) 


(2-10) 


°12  =  *  (A21  "  Al2)f 
or 

A12  =  b12  "  c12  » 
^1  =  °12  +  c12  • 

(B)  Ternary  Systems.  Wohl's  equation  of  the  third  order 

for  ternary  systems  is 

logtfi  ■   Z22  [A12   +  2Z1(A2lq1/q2   *  A3.2if 
«   Z32  [a13   +  2Z1(A3lq1/q3   -  A13)] 

•  Z2Z3  \h2\A\/<\2   ♦  AL3  -  A32q1/q3  »  2Z1(A3lq1/q3  -  Al3) 
+  2Z3(A32qi/q3  -  A23q1/q2)   -  C(l  -   2Z]  )J    .  (2-11) 

The  equations  for  logo2  and  log*JJ3  can  be  produced  by  a  rotation 

principle  from  the  subscript  triangle. 

A 

3* —  2 

In  order  to  obtain  the  equation  for  logl£2,  all  subscripts 
in  equation  (2-11)  are  rotated  one  unit  (  e.g.  x^-*x2,  A^2-*-A23) 
and  the  equation  for  log03  follows  from  the  equation  for  log#2  in 
the  same  way.  The  symbols  A^2t  A2^  and  the  effective  raolal  volume 
ratio  are  the  same  as  in  the  binary  equations.  The  ternary  cons- 
tants 0lt  C2  and  G3  which  have  to  be  determined  from  ternary  data 
are  related  by, 

G2  M  G1q2/q1    , 
G3  =  01q3/q1 

For  the  Margules  three-suffix  equation  G^=G2=G3,  since  qi/q2=l 

and  qi/q3=l. 

The  ternary  three-suffix  Margules  equations  can  be  obtained 


from  equation  (2-11): 

logKL  *   x22  [A12  ♦  **].(%  -  A12)J  +  x32  [al3+2x1(A3L-A13)] 
+::2x3  [a21   +  A13  -  A32  +  2xL(A31  -  A]3) 
+2x3(A32  -  A23)  -  G(l  -  2xL)J  .  (2-12) 

log02  and  log^3  can  be  obtained  by  the  rotation  principle.  These 
equations  a^e  identical  algebraically  with  the  three-suffix  eq- 
uations used  by  Benedict  et.  al.  (1)  and  ternary  Redlich-Kister 
equations  (13).  llala  (8)  has  shown  that  the  Redlich-Kister  equa- 
tion is  equivalent  to  the  Margules  equation  in  both  binary  and 
ternary  form. 

D^T^HMINATION  OF  CONSTANTS 

(A)  Binary  Marr<ules  Constants.  The  constants  A12,  A^  in 
Margules  equation  (2-7)  are  the  terminal  values  of  the  curves 
relating  the  logarithm  of  the  activity  coefficient  and  the  compo- 
sition of  the  solution: 

A12  =  lim  logtf,  , 

xi-»0  (3-1) 

A^  ■  lim  log#2  . 

x2->0 

The  values  of  the  constants  A12,  A2^  characterize  a  given  bi- 
nary system,  and  the  determination  of  those  constants  has  been  an 
important  problem  in  studies  of  phase  equilibria.  They  could  be 
determined  theoretically,  but  from  the  view  point  of  practical 
application  they  are  always  determined  from  experimental  data. 

The  calculation  of  the  constants  A-^  A^  from  the  known 
equilibrium  composition  of  the  vapor  and  liquid  is  the  most  precise 


method.  There  are  two  unknowns  A.2  and  Aj,  to  be  determined  in 
equations  (2-7).  Mathematically  speaking  equations  (2-7)  can  be 
considered  as  equations  in  the  two  unknowns  A^2  and  A21  which, 
therefore,  can  be  evaluated  for  each  data  point.   The  values  of 
A^2  and  A21  determine1  in  this  way  are  called  point  values  of  A12 
and  Ao^. 

Since  Margules  aquations  are  linear  in  the  constants, A^2  aa(^ 
Aji  they  can  also  be  determined  by  the  method  of  ordinary  least 
squares  or  the  weighted  least  squares  method  (  22,  23  ). 

In  this  work,  the  method  of  ordinary  least  squares  was  used  in 
evaluating  the  binary  Marseilles  constants.  This  was  due  to  the 
fact  that  no  information  was  available  regarding  the  precision  of 
experimental  measurements  and  therefore  it  was  assumed  that  all 
points  had  equal  wei;ht. 

The  constants  A^  and  A«^  were  calculated  by  computer  by  app- 
lying the  method  of  least  squares  to  four  different  arrangements 
of  equations  (2-7). 
Method  1. 

var.  1  =J[  loft(l2  -  A12  -  2(A21  -  A]_2)XlJ  2  ,  (3-2) 


^r^-xi 

Jl 


Var*   2  ^CtT^  -  A21  -  2<A12   -  A21)X2J  2      ^  (Jrt) 


Method  2. 
var 


.   1  =  Sjflog^  -    a-xx)2  ^2   -   2(A2X   -  A12)x1]J    2,      (3-<+) 
var.   2  =  y\}°^2  "   &**l>*  [^1  *  2(A12  "  A21^2fr    2,      (3-5) 


A.-I 


8 

Method  3. 

Var-=lLT^-Al2-2(A21-A12)Xl]2 

\£f0?£|fr  -  A2L  "  2(A12  -  A21)X2J2   f 

(3-6) 

Method  4. 

var.  =  Jl1o«^1  ~  <1-*1>2  C^  -  2(A21  -  A12)acj}2 

♦fflog^  -  d-x2)2  [A21  -  2(AL2  -  A21)x2]]2 

(3-7) 

Method  5. 

This  method  consisted  of  averaging  the  point  values. 

Method  1  is  equivalent  to  the  J-plot  used  by  Severns,  Sesonske, 

Perry  and  Pigford  (16),  or  to  Robinson  and  Gilliland's 

(14)  graph- 

ical  method.  Method  2  is  equivalent  to  plot  of  log#£ 

vs  X£.   For 

showing  the  differences  between  method  1  and  method  2, 

the 

following  substitutions  in  equations  (2-7)  were  made, 

al  =  A12  »                a2  =  A21  » 

bl  ■  2(A21  "  A12>  »        b2  =  2CA12  -  A2i)  , 
rx  =  log^  ,             r2  ■  log/)2  , 

(3-8) 

ax  m   (l-xi)2  ,           s2  =  (l-x2)2, 

it  follows 

rl  =  ^  al  +  blxl  ^81  • 

(3-9) 

r2  ■  (  a2  +  b2x2  )82. 

(3-10) 

Rearranging  eauation  (3-9),  (3-10) 

|^-  =»  ai  +  bixi  f 

(3-11) 

|2  =  a2  +  b2x2  , 
are  obtained. 

(3-12) 

Using  equations  (3-9)  and  (3-10)  to  obtain  a^,  b^,  a2,  b2 
by  the  method  of  least  squares  gives  values  different  from  those 
obtained  using  equations  (3-11)  and  (3-12).  This  can  be  shown 
by  the  following 

#1  •fiX*l  -   (a1   +  b1:c1)s])2   ,  (3-13) 

°K  s^&  -  (ai +  bixi^2 •  (3-w 

Minimum  values  of  &i   and  ^i   were  found  as  follows : 
By  partial  differentiation  with  respect  to  a^  and  b-L,  then  setting 
the  partial  derivatives  equal  to  zero,  the  following  equations 
were  obtained, 

^<*.  =  **\\*l  -   (ax  +  bix1)sj(-si)  ■  0   ,  (3-15) 

>bl   a  2£,Lr1  -   (ax  +  b1x1)sLj(-xi81)  =  0   ,  (3-16) 


__ 
2^ 


$i  =  *?£&  -   <aL  +  b1x1)](-l)  ■  0   |  (3-17) 


^,   ■  2£,[jg.  -   (a-,    +  blXl)J  (-xL)  =  0   ,                                      (3-18) 

Rearranging  equation   (3-15)  to   (3-18)  it  follows  that 

^-r^i    ■  ai^-si2   +  bi^xis-,2      ,  (3-19) 

^rjX^  =  a^x^  +  b^  x12s12    ,  (3-20) 

X|i  =  naL  +  bx^xj^     ,  (3-21) 

^*-BT  xl  =  al^xl   *  hX*!2     »  (3-22) 

Solving  equation   (3-19)  with  equation  (3-20)  ai  and  b^  were 
obtained, 


n  n 


al  ■      /^Iffl^ii2^    )  7iErlxT^(Sxls12)      ,  (3-23) 

*  (Isi*K  ix^s^;   -   (£x1siji)(2x1sL)  ♦ 

1  =      <nsl    )(^rTxlsl)   ~    (^x1s1)(^-x1s12)        .  (3-24) 

(gafeMgxffij*)  -   (^xls^X^xis^ 
a2*   ^2  can  bc  obtained  by  changing  subscripts  1  to  2, 

Also  solving  equation   (3-21)  with  equation   (3-22)   ai,   bi  were 
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obtained : 

al  -    ngx^  -  (gxj.X^) i  (3-25) 

»£fj-*L  -  (txi)(£|I) 
b*  =  n^V  -  (-Xl>(£Xli •  <3-26) 

a2  and  b2  were  obtained  by  chafing  subscripts  1  to  2, 

After  a^t  b]  ,  a2  and  b2  h4d  been  determined,  A^2  and  A21  w^*"® 
calculated,  but  A-^7  and  A,^  calculated  from  a^,  b^  were  not 
necessarily  equal  to  Aj2  and  Aj^  calculated  from  a2»  D2« 

Method  3  was  the  combination  of  equation  (3-2)  and  (3-3)  and 
method  4  was  the  combination  of  equation  (3-4)  and  (3-5). 

Methods  1  and  2  yield  two  sets  of  binary  constants  from  a 
binary  system  while  methods  3,  4  and  5  yield  only  one  set  of  binary 
constants.  A  detailed  description  of  the  computer  program  is 
presented  in  the  Append!  . 

(B)  Ternary  Constant  In  the  Margules  Equation.  Wohl  (27) 
defines  binary  constants  in  terms  of  interaction  constants 

\l  *   2al2  +  3a122 

*21  =  2a12  *  3aH2 

A23  =  2a23  *  3a233 

A32  *  2a23  +  3a223 

^3  "  2a13  +  3a133 

A31  =  2a13  +  3a113 
and  a  ternary  constant, 

C  =  3a112  +  3a133  ♦  3a223  -  6a123  ,  (3-28) 

where  a^^2»  a]_33»«»  are  ttxe  empirical  constants  which  measure  the 


(3-27) 
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interactions  in  various  group  of  molecules  112,  133... 

A  new  ternary  constant  C*  was  more  recently  defined  by  Wohl 
(28), 

C*  ■  7  (  a1L2  +  a122  +  au3  •  a133  +  a223  •  a233  )  -  6*l23- 

(3-29) 
By  apolication  of  equation  (3-27),  it  follows 

C*  -  C  ■  £(AL2  -  A21   +  A23  -  A32  +  A31  -  A13)   .       (3-30) 
In  this  work  the  older  definition  of  C  (27)  as  given  by 
equation  (3-28)  was  used. 

In  the  ternary  three-suffix  Margules  equation,  the  ternary 
constant  has  its  greatest  influence  on  the  calculated  value  of 
logOk  when  X£X.(l-2xjj)  has  its  maximum  value.  The  maximum  value 
of  xiXj(l-2xk)  occurs  for  very  small  values  of  x^  and  equal  values 
of  X£  and  xj;  the  term  vanishes  when  xj^  and  reverses  its  sign 
for  higher  values  of  xk.  The  most  suitable  region  of  concentratim 
in  which  to  determine  the  ternary  constant  from  experimental  data 
is  the  liquid  concentration  which  makes  the  value  of  X£Xj(l-2xk) 
as  large  as  possible.   In  this  work,  ternary  constants  were  deter- 
mined by  the  following  procedure. 

The  ternary  constant  G  was  calculated  from  the  binary  constants 
and  the  experimental  value  of  the  ternary  activity  coefficients 
from  the  following  equation, 

C  =  -  IoSV.x.  *  3-°g£al,  o-Q  tm  m  . 

xiXj(l-2xk)        •  (3-31) 

log"cal.  C»0  is  the  value  of  log J  calculated  with  only  binary 
constants.  Prom  a  given  set  of  experimental  data  three  different 
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average  values  of  C  were  evaluated  for  comparison: 

1)  using  all  data  points. 

2)  using  only  the  data  points  where  x  is  greater  than  0.05. 

3)  using  the  data  points  where  x  is  greater  than  0.05  as  well 
as  xiXj(l-2xjc)  greater  than  0.1. 

Concentrations  smaller  than  0.05  are  ignored  because  of  the 
difficulty  in  obtaining  accurate  analysis  of  the  component  when 
the  concentration  is  relatively  small  and  the  fact  that  suall 
errors  in  the  composition  reeult  in  large  variations  in  the 
activity  coefficient. 

The  average  values  of  the  ternary  constants,  which  were 
evaluated  under  the  condition  3,  were  used  to  obtain  a  more  satis- 
factory fit  in  the  activity  coefficient  equation. 

A  flow  diagram  and  a  detailed  computer  program  are  shown  in 
the  Appendix. 

RESULTS  AND  DISCUSSION 

A»  General  Results.   Vapor-liquid  equilibrium  data  for  only 
11  non-aqueous  ternary  systems  could  be  found  in  the  literature. 
Table  1  gives  the  components  of  these  systems  and  the  references. 

Binary  constants  used  in  the  calculations  were  obtained  by 
the  methods  described  in  a  previous  chapter  or  from  the  literature. 
These  are  listed  on  Table  2.  On  Figures  1  through  19  the  experi- 
mental data  and  the  Margules  equation  are  compared  for  most  of 
the  binary  pairs  comprising  the  eleven  ternary  systems.  On  these 
figures  the  solid  line  is  a  plot  of  the  Margules  aquation  with 
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constants  evaluated  by  method  4. 

esults  of  calculations  of  the  three-suffix  Margules  equation 
are  presented  in  Tables  4,  5,  6,  7,  3t  12,  13,  14,  15,  16  and  17 
which  give  the  ternary  constant  and  the  average  of  the  absolute 
value  of  the  percent  deviation  of  calculated^  with  C=0  from  the 
observed  Q  for  each  component  in  the  three  different  cases. 

1.  Using  all  data  points. 

2.  Using  only  the  data  points  where  x  is  greater  than  0.05. 

3.  Using  only  the  data  points  where  x  is  greater  than  0*05 
as  wall  as  x^xj(l-2xjt)  greater  than  0.1. 

When  it  was  found  necessary  to  use  a  ternary  constant  the 
absolute  value  of  the  percent  deviation  of  7f  calculated  with  C, 
from  the  observed  0  for  each  component  w|  also  tabulated. 

The  calculations  were  made  by  computer.   In  order  to  ensure 
that  no  mistakes  were  made  in  the  processing  or  programing,  the 
results  were  checked  by  hand  calculation  for  a  few  systems. 

B.  Discussion. 

System:  Benzene  -  Cyclohe:cane  -  Methylcello solve 

Figures  1  and  2  are  plots  of  logO  versus  x  for  the  Benzene  - 
Methylcellosolve  and  Cyclohexane  -  Methylcellosolve  binary  system*. 

The  three-suffix  Margules  equation  fit  the  data  in  Figure  1 
well  but  did  poorly  in  Figure  2.  The  absolute  value  of  error 
shown  in  Table  3  is  a  good  indication  of  how  well  the  data  are 
fitted  by  the  binary  Margules  equation. 


Table  1.   Components  of  Systems 


1>4 


System 
No. 

(1) 

Components  of 
(2) 

System 

(3) 

Ref. 

1 

Benzene 

'Jyclohexane 

tiiylcel  Losolve 

(20) 

2 

Benzene 

Gyclohexane 

Furfural 

(19) 

3 

Isooctene 

Toluene 

Furfural 

(21) 

k 

Benzene 

Gyclohexane 

Acetone 

(9) 

5 

MEK* 

n- Heptane 

Toluene 

(17) 

6 

Rfchgrlaoctata 

Benzene 

Gyclohexane 

(4) 

7 

Benzene 

Gyclohexane 

MIBK** 

(5) 

8 

Benzene 

Gyclohexane 

M^.K* 

(6) 

9 

Benzene 

Gyclohexane 

Trichloroetb/lene 

(12) 

10 

Benzene 

GClif 

Ethanol 

(2) 

11 

CCli, 

MEK* 

Cyclohexana 

(11) 

*  Methyl  ethyl  ketone 
**  Methyl  isobutyl  ketone 
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Table  2. 

Margules 

Binary 

Constants 

System 
No. 

A13 

A3I 

A3  3 

A32 

Al2 

A21 

Re- 
mark 

1 

0.3278 
0.3100 

0.7044 
0.6440 

0.7526 

0.6930 

1.1910 
0.9750 

0.1738 
0.1600 

0.1515 
0.1600 

a 
b 

2 

0.2450 
0.2419 

0.2830 
0.2528 

0.7420 
0.8812 

1.0460 
1.0260 

0.1600 
0.1738 

0.1600 
0,1515 

b 

a 

3 

0.9200 
0.884*: 

1.0137 
0.9214 

0.5232 

C.4832 

.    ?17 
0.3525 

0.1203 
0.1293 

0.1017 
0.1453 

a 
b 

4 

0.1107 

0.1773 

0.8420 

0,4010 

0.1580 

0.1580 

b 

5 

0.1995 
0.1942 

0.1184 
0.1202 

0.1422 
0.1430 

0.1047 

0.U 

0.4965 
0.4894 

0.5150 
0.470*-. 

a 
b 

6 

0.4253 

0.3621 

0.1364 

O.i: 

0.041C 

C.0328 

a 

7 

0 
0 

0 
0 

0.1000 
0.0690 

0.4100 
0.4421 

0.158G 
0.1364 

3.15S0 
0.1567 

b 

a 

8 

0.1204 

0.2685 

0.468' 

0.4546 

0.1364 

0.1567 

a 

9 

0 

0 

0 

0 

0.1100 
0.1026 

0.1800 

0.168--. 

1.1530 
0.1364 

0.1580 
0.1 

b 

a 

10 

0.5850 

0.8510 

0.6240 

0.9570 

0.0420 

0.0337 

a 

11 

0 

0 

0.5010 

0.4040 

0.1 600 

0.2200 

b 

a 

Calculated 

by  method 

4. 

b 

Prom  liters 

ture. 

16 


0.5  1.0 

MOLE     FRACTION   BENZENE    IN    LIQUID 
FIG.  I    BENZENE-METHYLCELLOSOLVE    SYSTEM 
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CYCLOHEXANE 
METHYLCELLOSOLVE 
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MOLE    FRACTION     CYCLOHEXANE     IN    LIQUID 
FIG.2    CYCLOHEXANE-METHYLCELLOSOLVE    SYSTEM 
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Table  3,  Average  of  the  Absolute  Value  of  the  Percent 

Deviation  of  (fcax  from  tf^xp  for  binary  Systems 


Binary 

A12 

A21 

Absolute  value  of 
percent  error 

systems 

Component  1 

Component  2 

Benzene  - 
Methylcellosolve 

0.3273 

0.7044 

2.8% 

4.5% 

Cyclohexane  - 
Methylcellosolve 

0.7526 

1.1910 

10. 

6.4% 

Benzene  - 
Cyclohexane 

0.1738 

0.1515 

0.2% 

0.3% 

Table  4.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  Ternary  ocal  from  o&xp 


Component 

Binary  constants 
from  Reference  21 

Binary  constants  by 
method  4 

C  =  0 

C  =  -1.6 

C  =  0 

C  =  -1.5 

Benzene 

Cyclohexane 

Methylcellosolve 

15.96 
16.31 
13.43 

6.85 

13.21 

6.28 

15.50 

12.60 

8.18 

5.45 
13.03 
10.37 

A  comparison  of  the  absolute  value  of  the  deviation  of  *cai 
with  and  without  C  from  0exp  is  shown  in  Table  4.   It  appears 
that  the  system  needs  a  ternary  constant  for  more  satisfactory 
ternary  activity  coefficient  calculation.  This  might  be  due  to 
the  fact  that  only  a  few  experimental  points  were  presented  for 
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the  binary  system  (3yclohe:cane  -  Methylcellosolve  and  the  Margmles 
binary  equation  did  not  correlate  the  binary  data  very  well.  On 
the  other  hand  the  other  two  binary  systems  were  fitted  quite 
well  with  the  Margulcs  equation*  The  agreement  is  better  using 
Thornton  and  Garner* s  constant!  than  with  computer  calculated 
constants.  Also  their  ter.Tiary  constant  made  the  calculated  results 
more  satisfactory. 

System:  Benzene  -  Tyclohcxane  -  Furfural 

Plots  of  logO  versus  x  for  the  three  binary  systems,  Benzene  - 
Cyclohexane,  Benzene  -  Furfural,  and  Tyclohexane  -  Furfural  are 
shown  in  Figures  3t  bt   and  5,  The  binary  data  for  these  three 
binary  systems  fit  the  binary  Margule  I  ■      ns  fairly  well* 

Table  5  shows  the  absolute  value  c      percent  error  of  the 
calculated  activity  coef     its  from  the  experimentally  deter- 
ged activity  coefficients,  Thes-^     eonpared  With  and  with- 
out  use  of  a  ternary  constant.  The  ternary  constant  does  not 
improve  the  results  in  activity  coefficient  calculation  very  much. 
It  is  evident  that  a  ternary  constant  is  not  necessary  in  the 
system. 
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Tabic  5.  Average  of  the  Absolute^ Value  of  the  Percent 
Deviation  of)  cal  from  0  ^-p 


Component 

Binary  constants 
from   Reference  21 

Binary  constants  by 
method  4 

G  =  0 

3=-0.251 

C  =  0 

G=-0.251 

Benzene 

8.14 

.17 

9.60 

7.14 

Cycloho^nc 

7.17 

8.93 

9.08 

12.11 

Furfural 

5.15 

4.65 

7.33 

7.21 

System:   Isooctane  -  Toluene  -  Furfural 

Figures  6t  7  and  8  are  plots  of  logo  versus  x  tor  Isooctane  - 
Toluene,  Toluene  -  Furfural,  and  Isooctane  -  Furfural.  The  binary 
Margules  equation  correlates  the  binary  systems  reasonably  well. 

Table  6  shows  t:     e  ternary  constant  C=-0.317  proposed 
by  Thornton  and  Garner  is  not  much  help  in  correlation  of  the 
activity  coefficients  of  the  system.  The  small  value  of  G  in- 
dicates that  activity  coefficients  in  the  ternary  system  could  be 
correlated  relatively  well  without  the  ternary  constant. 


Table  6.  Average  of  the  Absolute /Value  of  the  Percent 
Deviation  of  & cai  from  lexp 


Binary 

constants 

Binary 

constant  by 

>omDonent 

from  'teferenoe  21 

method  4 

G  =  0 

G=-0.317 

G  =  0 

C=-0.200 

8.17 

4.82 

G.99 

5.09 

Isooctane 

3.38 

3.19 

5.76 

4.18 

Toluene 

3.50 

3.29 

4.36 

4.63 

Furfural 
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MOLE    FRACTION     ISOOCTANE     IN     LIQUID 
FIQ.6     ISOOCTANE-TOLUENE      SYSTEM 
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MOLE     FRACTION     ISOOCTANE    IN    LIQUID 
FIG.  8     ISOOCTANE- FURFURAL    SYSTEM 
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System:  Benzene  -  (Jyclohexane  -  Acetone 

Binary  constants  were  not  calculated  for  this  system  because 
the  binary  -quilibriura  data  were  reported  only  in  the  form  of 
constants  in  the  three-suffix  Margulos  equation* 

The  to.mr.ry  constant,  G*",  which  the  investigators  of  this 
system  used  in  c  uivalent  to  A21  +  A13  +  A32  -  G.  (  i.  e. 
G  =  A21  ♦  A^3  +  A32  -  Cfc  ).  Using  their  reported  value  of  C  and 
the  above  relation,  a  value  of  -0.03  was  obtained  for  C.  Values 
of  w1#  °2»  "3  were  calculated  with  and  without  this  ternary  cons- 
tant. The  deviation  of  these  calculated'  froin  experimental  7 
are  saown  on  Table  7*   It  can  be  se^n  that  a  ternary  constant  is 
not  necessary  in  the  three- suffix  Margules  eq     is.  Further- 
more  Table  7  shows  that  if  the  dot  points  which  satisfy  ttie 
condition  x  smaller  than  0*05  were  omitted,  the  absolute  value  of 
the  average  deviation  would  be  reduced * 

Table  7.  Av     of  the  Absolute jValue  of  the  Percent 
Deviation  of  »cai  from  X^m 


Component 

All  data  points 

The  points 
x>J 

where 
.05 

No.  of  dat     Deviation 

No.  of  data 
point ~ 

Dev. 

points    0  ■  9 

C=-0. 03 

C  =  0' 

Benzene 

Cyclohcrrane 

Acetone 

'0       7.47 
60      14.76 
60        6.63 

7.99 

14.20 

6.69 

52 
46 
45 

.JO 

11.50 

5.65 
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Systems:  Methyl  ethyl  ketone  -  Heptane  -  Toluene 
Table  8  shows  that  the  average  of  the  absolute  value  of  the 
percent  error  of'  for  X3  greater  than  0.05  is  very  much  reduced 
from  the  average  absolute  percent  error  obtained  from  all  data 
points.  This  is  probably  due  to  the  difficulty  in  obtaining 
accurate  analysis  of  the  component  present  in  relatively  small 
concentrations.  Figures  9,  10,  and  12  show  that  the  three  cons- 
tituent binaries  of  the  ternary  system  are  excellently  represented 
by  the  binary  Margules  equation.  The  average  of  the  absolute 
value  of  the  percent  error  between  calculated  and  experimental  " 
for  these  three  binary  systems  is  less  than  2  percent.  These 
systems  are  listed  on  Table  9,  The  Margules  constants  .eterrainei 
by  several  methods  are  very  close  to  each  other.  They  are  listed 
on  Table  10. 

Table  8.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  "o  for  Ternary  System 


Binary  constants 
from  Reference  17 

Binary'  constants  by 
method  4 

Component 

All 

data 

:c>0.05 

All 

\:a 

x>0.05 

9?* 
data 

dev. 

of   dev, 
data 

data 

dev. 

of   dev. 
;ata 

Methyl- 

ethyl- 

ketone 

60 

2,10 

52   1.56 

60 

2,76 

52   2.22 

Heptane 

60 

2,20 

46   1.70 

60 

2.55 

46   2.12 

Toluene 

60 

14,50 

45   1.34 

60 

15.75 

45   2.12 

The  ternary  system  shows  less  than  3  percent  deviation  be- 
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tween  calculated  b   and  observed  "  t  therefore  no  ternary  constants 
are  necessary. 


Table  9,  Absolute  Value  of  the  Percent  Deviation 
of  *cal  from*exp 


MEK*Toluene 

n-Heptane-Toluene 

M^K-Ifeptane 

A13    A31 
0.1995  0.1184 

A23    A32 
0.1^22  0.1047 

A12     A21 

0.4965   0.5150 

Component  1 

0.56 

0.73 

1.26 

Component  2 

0.37 

0.62 

1.21 

Table  10.  Margules  Binary  Constants  for  the  System 
MEK  -  n-Heptane  -  Toluene 


A13 

A31 

A23 

A32 

A12 

A21 

Method  1 

0.1809 
0.2100 

0.1284 
0.0941 

0.1207 
0.1710 

0.1166 
0.0811 

0.3799 
0.5286 

0.5293 
0.4240 

a 
b 

Method  2 

0.1811 
0.1990 

0.1277 
0.1168 

0.1381 
0.1^64 

0.1114 
0.9993 

0.4658 
0.5181 

0.5053 
0.4557 

a 
b 

Method  3 

0.2104 

0.1280 

0.1711 

0.1154 

0.5214 

0.5371 

Method  4 

0.1995 

0.1184 

0.1422 

0.1047 

0.4965 

0.5150 

Method  5 

0.1983 

0.1172 

0.1-'24 

0.1029 

0.4927 

0.5053 

From 
Reference 

17  0.1942 

0*1202 

0.1 '+30 

0.0950 

0.4894 

0.4706 

a    determined  from  |"x  data  for  component  1. 
b    determined  from  V*  data  for  component  2. 


System:  ^thylacctate  -  Benzene  -  Cyclohexane 

The  binary  constants  determined  by  method  4  describad  in  the 
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previous  section  are  in  agreement  with  those  of  Qhac  and  Hougen 
(4)  leterroined  by  the  method  of  least  squares  using  a  computer. 
Their  tr.To  constant  Tledlich-Kister  equation  ifl  equivalent  to  the 
three-suffi::  Ma-gulcs  equation  of  Wohl.  The  constants  are  present- 
ed on  Table  11  for  comparison. 

Table  11.  Comparison  of  Binary  Constants 

From  Reference  13 
Binary  System    R,*K  Constants      Ilargules   llargules  Constants 

,  Constants      by  method  4 

b!2      c!2    A12     A21     A12     A21 

Ethylacetate(l)  0.03r>9  -0.0041  0.0410  0.0328  0.0410  0.0328 
Benzene  (2) 

Bcnsene  (1)      0.1443  0.0100  0.1343  0.1543  0.12G4  0.1558 
Cyclohexane (2 ) 

Cyclohexane(l)   0.3849  -n.0321  0.4170  0.3528  0.4253  0.3r>21 
Hthylacetate(2) 

Figures  11,  13  and  14  are  plots  of  log  ft  versus  composition 
for  the  three  binary  systems.  Agreement  is  very  good  over  all 
concentrations.   The  ternary  three- suffix  llargules  equation 
without  the  ternary  constant  gives  very  good  results  in  correlation 
of  ternary  vapor-liquid  equilibrium.  The  results  are  shown  in 
Table  12. 
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Table  12.  Average  of  the  Absolute/  Value  of  the  Percent 
Deviation  of  6cal  from  oexp 


Component 


ah    daffl  the  data  where     the  date  where 

RLL  aata  x>0.05  x>0.05,  u<£.l 


Ethylacetate  77      1.14      77  1.14      14     1.37 

Benzene      77      0.87      76      0.323     11     1.37 
Cyclohexane  77  0.917     77      0.917     17     1.75 


This  is  also  in  agreement  with  Chao  and  Hougen's  statement, 
(4)  "No  additional  constants  were  require!  for  ternary  effects.'* 

System:  Benzene  -  Cyclohexane  -  Mathyliscbutylketone 

For  the  Benzene  -  Cyclohexane  system  the  binary  constants 
were  taken  from  the  previous  work.   Benzene  -  MetLylisobutylketone 
is  nearly  an  ideal  solution  and  xero  can  be  used  for  these  binary 
constants. 

The  binary  constants  for  the  Cyclohextne-Methylisobutylketone 
system  were  determined.  The  plot  of  log  •  versus  composition  is 
given  in  Figure  15,  An  exanination  of  the  Figure  15  reveals  that 
this  data  did  not  fit  the  curve  well  when  x^  approaches  zero, 
this  poor  correlation  might  cause  more  deviation  when  the  binary 
constants  of  this  system  are  used  in  the  ternary  three-suffix 
Margules  equation  to  predict  the  vapor-liqxiid  equilibrium. 
Deviations  of  the  activity  coefficients  in  the  ternary  system 
are  presented  in  Table  13. 
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Table  13.  Average  of  the  Absolute, Value  of  the  Percent 
Deviation  of  'ca^  from  o  ( 


exp 


Component 

Binary  Constants 
from  Reference  5 

Binary  Constants  by 
method  4 

c  =  6 

C=-0.26 

C  ■  0 

C=-0.118 

Benzene 

3.56 

3.98 

3.59 

3.48 

Cyclohexane 

9.52 

8.37 

9.71 

8.82 

MIBK 

6.04 

7.47 

6.17 

6.50 

In  reference  to  Table  13,  the  ternary  constant  reduced  the 
deviation  of  component  two,  but  increased  the  deviation  of  the 
other  two  components.   From  the  view  point  of  the  total  devia- 
tion of  three  components,  the  ternary  constant  is  not  necessary 
in  this  system.   It  seems  that  no  ternary  effect  exists.  The 
above  discussion  is  not  in  agreement  with  the  statement  of  Rao 
et.  al  (5).  They  state  "  As  a  first  approximation  C  was  assumed 
to  be  zero  and  the  activity  coefficients  were  calculated;  good 
agreement  was  not  obtained  between  the  calculated  and  experimental 
values.  In  order  to  obtain  a  more  satisfactory  fit  of  the  data.. 
C  was  calculated  ....  By  inserting  these  seven  constants  into  the 

three-suffix  Margules  equation  of  Wohl The  average  error 

between  calculated  value  and  experimental  activity  coefficient 
in  ^.  ,  02  anc^  ^3  arc  respectively  2.4  percent,  4  percent,  and 
4.5  percent.  " 

This  disagreement  might  be  lue  to  the  fact  that  the 
calculations  performed  at  K.S.U.  were  done  by  computer,  whereas 
their  calculations  were  probably  done  by  hand.  In  this  work 
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average  deviations  are  reported  as  absolute  values;  deviations 
reported  by  Rao  et.  al.  (5)  were  most  likely  averaged  algebraically. 
System:  Benzene  -  Cyclohexane  -  Metuylethylketone 
These  three  binary  systems  were  correlated  very  well  by  the 
three-suffix  Margules  equation.  Figures  16  and  17  show  plots  of 
log^  against  x  for  two  of  the  binary  systems.  The  plot  for 
Benzene  -  Cyclohexane  has  already  been  presented  in  Figure  11. 
In  Table  14,  the  deviations  for  this  system  are  seen  to  be  within 
limits  of  experimental  accuracy.  Hence  the  three-suffix  Margules 
equations  represent  the  data  fairly  well  without  the  ternary 
constant. 

Table  14.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  ^ca]_  from  tfexp 


Component 

All  data 

The  data  where 
x>0.05 

The  data  where 
x>0.05,  u>0.1 

No.  of 
data 

Dev. 

NSataf    Dev- 

NSataf    **■ 

Benzene 

86 

5.47 

77      4.62 

5     2.31 

Cyclohexane 

86 

6.11 

78      5.77 

14     9.51 

MEK 

86 

7.98 

81      7.16 

13    11.75 

System:  Benzene  -  Cyclohexane  -  Trichloroethylene 

Table  15  shows  clearly  that  for  this  system  the  small  ternary 

constant  was  of  no  use  in  the  Margules  three-suffix  equations. 

Furthermore  the  binary  constants  exerted  more  influence  than 

the  ternary  constant  in  predicting  vapor-liquid  equilibrium  in 

this  ternary  system. 
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FIG.  16      BENZENE-MEK     SYSTEM 


1.0 


37 


j_ 


X 


_L 


_L 


_L 


-L 


JL 


_L 


0.5  1.0 

MOLE     FRACTION     CYCLOHEXANE  IN    LIQ. 
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In  the  determination  of  binary  constants,  the  Benzene  - 
Trichloroethylene  system  behaves  ideally,  hence  these  binary 
constants  are  zero. 

The  constants  of  the  well  known  system,  Benzene  -  Gyclohexane 
were  previously  determined.  For  the  Gyclohexane  -  Trichloro- 
ethylene system,  log  0  versus  x  is  plotted  in  Figure  19.  The 
binary  Margules  equation  fits  the  experimental  data  fairly  well. 

Table  15.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  0  cal  from  l(exp 


Component 

Binary 
from  Re 

C  =  0 

Constants 
ference  11 

C=0.094 

Binary  Constants  by 
method  4 

C  ■  0 

00.094 

Benzene 

2.60 

2.69 

2,64 

2.69 

Gyclohexane 

3.20 

3.28 

2.93 

3.35 

TGE 

6.38 

5.73 

6.24 

5.50 

System:  Benzene  -  Carbon  tetrachloride  -  Ethyl  alcohol 
Figures  18,  20  and  21  show  plots  of  logi   versus  composition 
for  the  three  constituent  binary  systems.   It  was  found  that  the 
three-suffix  Margules  equation  correlated  the  three  binary  sys- 
tems well  except  for  a  few  points  which  were  obviously  due  to 
experimental  errors. 

Using  six  binary  constants  in  the  three-suffix  Margules 
equation,  the  activity  coefficients  were  calculated  for  the 
Benzene  -  Carbon  tetrachloride  -  Ethyl  alcohol  system  and  the 
deviations  from  the  observed  activity  coefficient  are  listed 
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FIG.   18      ETHANOL-BENZENE      SYSTEM 
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in  Table  16. 

It  appears  that  no  ternary  constant  is  necessary  in  this 
ternary  system.  Excluding  concentrations  less  than  0.05  the 
agreement  between  calculated  0  and  observed  6  was  found  even 
more  satisfactory. 

Table  16.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  o  cal  f rom  •  exp 


Component 

All 

data 

The  data  where 
x>0.05 

The  data 
x>0.05. 

No.   of 
data 

where 
u^O.l 

No.   of 
data 

Dev. 

No.   of        ^ 
data         LH5V# 

>ev. 

Benzene 

51 

8.99 

44            6.16 

37 

4.33 

CCl^ 

51 

5.93 

36            4.55 

31 

3.56 

Ethyl  alcohol 

51 

8.67 

41            4.00 

39 

4.07 

System:   Carbon  tetrachloride-Methylethylketone-Cyclohexane 
No  binary  data  have  been  reported  for  these  systems.   How- 
ever binary  Margules  constants  are  available.  The  ternary 
activity  coefficients  were  computed  by  substituting  the  six 
binary  constants  and  one  ternary  constant  given  by  Rao  (11) 
into  the  three-suffix  Margules  equation.  The  leviation  from  the 
observed  activity  coefficients  are  listed  in  Table  17. 
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Table  17.  Average  of  the  Absolute  Value  of  the  Percent 
Deviation  of  i ca]_  from  ^exp 


Component 

Deviation 
computed  by  author 

Deviation 
reported  in  reference 
11 

C  =  0 

C=-0.314 

C=-0.314 

CCl^ 

5.99 

3.28 

3.54 

IT]K 

10.14 

8.03 

4.43 

Cyclohexane 

6.62 

5.83 

4.24 

In  this  work  the  absolute  value  of  the  deviations  were  taken; 
whereas  Rao  did  not  mention  whether  his  deviations  were  absolute 
or  not.  In  Table  17  all  calculations  were  male  with  Rao's  report- 
ed binary  constants.  The  system  seems  to  need  a  small  ternary 
constant  to  make  the  ternary  Margules  equation  fit  the  data  more 
satisfactorily. 
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CONCLUSION 

The  results  show  that  if  each  of  the  three  binary  systems 
comprising  the  ternary  system  could  be  correlated  well  by  the 
three-suffix  Margules  equation,  the  prediction  of  the  ternary 
relationships  could  be  made  without  a  ternary  constant  in  most 
cases.  The  precision  of  the  prediction  of  an  activity  coefficient 
in  the  ternary  system  is  about  the  same  as  in  the  binary  Margules 
equation  fit  of  the  binary  data.  For  example,  in  the  Ethyl- 
acetate  -  Benzene  -  Cyclohexane  system  each  binary  Margules  eq- 
uation fit  the  data  to  within  one  percent,  and  the  ternary  activity 
coefficients  were  pre  'icted  to  about  one  percent  error.  This 
should  not  be  the  case  in  a  ternary  system  where  one  of  the  three 
constituent  binary  systems  was  not  well  represented  by  the 
Margules  equation. 

The  activity  coefficient  is  sensitive  to  temperature,  but  the 
three-suffix  Margules  equation  still  can  be  applied  to  isobaric 
systems  where  the  temperature  range  in  the  ternary  systems  is 
approximately  the  same  range  as  in  the  binary  systems.  The  Methyl - 
ethylketone  -  n-Ifeptane  -  Toluene  system  is  a  good  example  of  this. 

Five  different  ways  of  determining  Margules  binary  constants 
were  presented.  As  far  as  the  computer  calculation  is  concerned, 
method  4  is  the  best  one.  The  results  also  show  that  method  4 
is  more  satisfactory  in  ternary  activity  coefficient  calculation. 

The  results  indicate  that  systems  which  were  fitted  well 
without  a  ternary  constant  usually  showed  small  deviations  from 
ideal  behavior. 
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NOMENCLATURE 

A 

Binary  constant  (  double  subscripts  );  Antonine  constant 
(  single  subscripts  ) 

a 

Empirical  constant 

B 

Antonine  constant 

b 

Redlich-Kister  constant 

C 

Ternary  constant;  Antonine  constant 

c 

Redlich-Kister  constant 

G 

Total  free  energy 

n 

Number  of  data  points;  Number  of  moles  of  component 

P 

Vapor  pressure 

q 

Effective  molal  volume 

R 

Gas  constant 

r 

The  function  defined  in  Equation  (3-8) 

s 

The  function  "efined  in  Equation  (3-8) 

T 

Absolute  temperature,  °K 

3C 

mole  fraction 

Z 

Effective  molal  volume  fraction 

A 

Increment 

AGE 

Molar  excess  free  energy 

In 

logarithm  base  e 

log 

logarithm  base  10 

var 

Sum  of  squares  of  deviations  of  measured  values  from 
calculated  values 

i 

Activity  coefficient,  standard  state  taken  as  pure  sub- 
stance 
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9-  The  function  lefined  in  Equation  (3-13) 

(jy  The  function  defined  in  Equation  (3-14) 

Subscripts 

1  Component  1 

2  Component  2 

3  Component  3 

i  i**"  component 
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APPENDIX 


Description  of  Computer  Program 


a) 

Title:   Program  for  the  calculation  of  Margules  constants  by  the 
method  of  least  squares. 

Computer:   IBM  1620  with  60,000  positions  of  storage. 
Program  language:   Fortran. 

Description:  The  program  was  used  to  determine  binary  Margules 
constants  by  the  method  of  least  squares.  A  maximum  of  50  data 
po5.nts  was  provided.  The  point  values  of  A]_2  and  A21  at  each 
liquid  composition  were  calculated  first.  The  binary  constants 
were  calculated  by  five  methods  then  these  constants  were  used  in 
equation  (2-7)  and  the  average  of  the  absolute  value  of  the 
percent  deviation  between  calculated  and  experimental  activity 
coefficient  was  evaluated  for  each  of  the  five  methods. 

Input  cards  were  read  in  the  following  order; 

Card  1.  Title  or  comment  less  than  49  characters  including 
blanks  in  the  body  of  the  comment. 

Card  2.  M,  N 

M  is  a  fixed  point  constant  showing  the  type  of  input  activity 
coefficient  data.  N  is  also  a  fixed  point  constant  representing 
the  number  of  sets  of  experimental  measurements. 

Sample  input  data  cards  were  arranged  as  following  depending 
on  a  fixed  point  constant  M; 
for  M  =  1   Constant  pressure  with  x,  y,  t 
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Card  3,  A^t  Blt  Glf  A21  B2,  G2 

Card  4.  IT 

Card  5.  xlf  ylf  ti 

Card(4+N)xN,  yN,  tN 
for  M  »  2    constant  pressure  with  x,  y,  pi,  P2 
Card  3.     TT 
Card  4.     xlt  ylt  pu,  P21 

Card  (3+N)  Xj,,  yN,  p1N,  p2N 
for  Mi  3     constant  temperature  with  x,  y,  TT 
Card  3.     p^0,  p2° 
Card  4.     xj_,  yx,  ^1 

Card  (3+N)  xN  yN,  fl"N 
for  M  ■  4     t    are  available 

Card  3.     xlt  ylt  in,  ^21 


• 


Card   (2+N)     xN,  yN,  l(lN>  #2N 
for  M  =  5         In  &     ai*e  available 

Card  3.  xlt  y^t    (ln^)^,    (ln«2>l 

Card  (2+N)  Xjj,  yN,  (ln?i)Nf  (ln02)N 
for  M  =  6    logl  are  available 

Card  3.     xlf  ylf  (logYj.)^  (Iog0  2)l 
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Carl  4.     x2,  y2»  (logii>2»  (log*2>2 


Card  <2+N)  xN,  yN,  (lo^6i)N,  (logYo)^ 
All  output  is  printed  on  the  console  typewriter  in  the 
following  format: 

Title:   Benzene  -  Gyclohexane 
The  point  value  of  Margules  constants, 
xL        A12        A2i 


Average 
Benzene 


Method  1 
Method  2 
Method  3 
Method  4 
Method  5 

Gyclohexane 


A12 


A21 


Percent  Error 


A12 


^1 


Percent  Error 


Method  1  

Method  2  

Method  3  

Method  k  t 

Method  5  

Console  setting  for  operation: 

Overflow  check  switch  -  program 

All  other  check  switches  -  stop 

Sense  -witches  1,  2,  3,  k       immaterial 
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BLOCK     DIAGRAM     FOR    USE    WITH     IBM     1620 
FOR    COMPUTING     MARGULES    CONSTANT 


START 


I 


READ.  M  ,  N 


READ,AI  ,BI 
Ct,A2,B2fC2 


READJT 


±3. 


READ,££ 


READJT 


X    Y   T 


CALC. 
7»  log? 


± 


I 


READ, 
X,Y,TT 


X  ,  Y 

P     P 


CALC. 


CALC. 
"/»  og"" 


READ 
X,Y 

1,Sl 


CALC. 
log  7 


CALC. 
T,logT 


READ 
X,Y 

log7»ilogYa 


CALC. 

y 


• 

9 

S3 

CALCULATE     POINT 

VALUES    OF  A,2,A2| 

i 

1 

CALCULATE    Ai2,Azi 

BY    EQ.O-2)     TO 

. 

EQ.(3-7) 

» 

CALCULA1 

i 

E    AVERAGE 

VALUE    OF   ABSOLUTE 

DEVIATION     USING 

DIFFERENT    AI2.A2I 

1 

PRINT  ,  Al2,A2l 

PERCENT    ERROR 

FOR     5     METHODS 

• 

1 

• 

START 

s 
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G  Dr      'ATION  OF  MARGULES  CONSTANTS  BY  THE  METHOD  OF  LEAST  SQUARES 
DIMENSION  SXU(2),  GAL0C(2,50),  U(2,50),  SN(2)t  SU(2),  SX(2) 
DIMENSION  BA1(2,50),  BA2(2,50),  BA(2,50),  SBA(2),  AB(2),  SW(2) 
DIMENSION  SVR(2\  SRR(2),  SUV(2),  SUR(2),  GAA(2f50),  GAB(2,50) 
DIMENSION  SWX(2),  SWXX(2)g  S  V(2),   GVX(2),  0C(2)f  CD(2)t  GAC(2,50) 
DIMENSION  VAPRO(2),  CE(2),  CF(2),  GCLOC(2,50) ,  GDLOC(2,50) 
DIMENSION  A(2),  B(2),  C(2),PO(2) ,X(2,50) ,Y( 2, 50),T(50) ,S1(2) 
DIMENSION  GA(2f50)tGALN(2,50),GALOE(3,50),SXa:(2),AA(2),B3(2) 
DIMENSION  S2(2),A12(2),A21(2),DEV(2,50),P(2,50),S4(2),S5(2) 
DIMENSION  GBLOC(2,50),S3(2),SM(2),GAD(2,50),SV(2),GAE(2,50) 
DIMENSION  GEL0C(2,50) 

1  FORMAT  (E12.6) 

2  FORMAT  (12,  12) 

3  FORMAT  (  E12.6,  E12.6  ) 

4  FORMAT  (  EI2.6,  E12.6,  73.2.6,  E12.6  ) 

5  FORMAT  (  E12.6,  E12.6,  E12.6,  E12.6,  E12.6  ) 

6  FORMAT  (  E12.6,  E12.6,  E12.6,  E12.6,  E12.6,  E12.6  ) 

7  FOR7!AT(6HMETHOD  14,  2X  E12.6,  E12.6,  EJL2.6  ) 

8  FORMAT  (12HAVERAGSVALUE  El  .  ,   I  .   ) 

9  FORMAT  (A7X  3HA12  9X  3HA21  6X  10H     NT  ER  ) 

10  PAUSE 
READ  16 
REA^  2,  M,  N 
PRINT  16 

GO  TO  (  11,  12,  13,  14,  14,  14  ),  M 

11  READ  6,  A(l),  B(l),  C(l),  A(2),  B(2),  G(2) 

12  READ  1,  PI 
TO  TO  14 

13  3,  VAPRO(l),  VAPR0(2) 

14  SO--0. 

so=o. 

SRsO. 

SS=0. 

SW=0. 

DO  15  1=1,2 

S1(I)=0. 

S2(I)=0. 

S3tt)=0. 

S4(I)=0. 

S5(I)=0. 

33A<I)=0. 
SGV(I)=0. 
SGVX(i)=0. 
SN(I)=N 
SM(I)=N 
SU(I)=0. 
SV(I)=0. 
Cl)=0. 
SXU(I)=0. 
SRR(I)=0. 
SVR(I)«0- 


1 

53 

SW(I)=0. 

swx(i)=o. 

swxx(i)=o. 

SUR(I)=0. 

SUV(I)=0. 

15 

SXX(I)=0. 

16 

FORMAT  (49H 

) 

17 

FORMAT  (38HTHF.  POINT  VALUE  OF  MARGUL3S  CONSTANTS) 

18 

MAT  (/5X  2HX1  9X  3HA12  9X  3HA21 ) 

19 

FORMAT  (  E12.6,  E12.6,  E12.6  ) 

DO  20  J=1,N 

READ  5,  X(1,J),  Y(1,J),  T(J),  P(1,J),  P(2,J) 

Y(2,J)=1.-Y(1,J) 

20 

X(2,J)=1.-X(1,J) 

DO  33  1=1,2 

DO  33  J=1,N 

GO  TO  (  21,  23,  22,  24,  26,  28  ),  M 

21 

P(I,J)=EX?(2.302535*(A(I)-B(I)/(C(I)+T(J)))) 
GO  TO  23 

22 

P(I,J)=VAPRO(I) 
PI=T(J) 

23 

GA(I,J)=PI*Y(I,J)/P(I,J)/X(I,J) 
GO  TO  25 

24 

GA(I,J)=P(I,J) 

25 

GALN(I , J)=L0G(GA(  7. ,  J)  ) 
GO  TO  27 

26 

GALN(I,J)=P(I,J) 
GA(I,J)=EXP(GALN(I,J)) 

27 

GALOE ( I , J )=GALN ( I , J )/2 . 30258 5 
GO  TO  29 

28 

GALOE(I,J)=P(I,J) 

GA(I , J)aEXP(2. 302585*GALOE(I , J) ) 

29 

U(I,J)=GAL0E(I,J)/((1.-X(I,JJ**2) 

30 

FORMAT  (///) 

32 

SU(I)=SU(I)+U(I,J) 
SX(I)»SX(I)*X(ItJ) 

SXX(I)=SXX(I)+X(I,J)*X(I,J) 

sxu 

SV(I)=SV(I)+(1.- C(I,J))**4 

SW(I)=SW(I)  +  (1.-2.*X(I,J))**2 

SVR(I)=SVR(I)+(1.-2.*X(I,J))*2.*X(I,J) 

3RR(I)=3RR(I)+'+.*v(I,J)**2 

SUV(I)aSUV(I)+U(IfJ)*(l,-2.*X(I,J)) 

S:JR(I)=3UR(I)+U(I,J)*2.*.(I,J) 

SWX(I)sSWX(I)  +  (l.-X(I,J))**4*X(I,J) 

swn(i)=swxx(i)+(i.-x(i,j))**if*x(i,j)**2 

SGV(I)=SGV(I)+GALOE(I,J)*(l.-X(I,J))**2 

SGVX(I)sSGVX(I)-»-GALOE(I,J)*(l.-X(I,J))**2*;v(I,J) 

33 

CONTINUE 
PRINT  30 
PRINT  17 
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DO  50  1=1,2 

DO  49  J=1,M 

GAL0"(3,J)=GAL0R(1,J) 

BAl(I,J)aGALOR(I,J)/(l.-X(I,J))**2 

BA2(l,J)«2.*X(I,J)*(GALnE(I+l,J)/X(I,J)**2-BAl(I,J)) 

BA(I,J)=BA1(I,J)+BA2(I,J) 

IF   (   X(I,J)-0.i;   )   35,    36,    34 

34  IF   (   X(I,J)-0.Q5    )   36,    36,    35 

35  SN(I)=SN(l)-l. 
GO  TO  49 

36  SBA(I)=SBA(I)+BA(I,J) 

49  CONTINUE 

50  SBA(I)=S3A(I)/SN(I) 
PRINT  18 

DO  60  J=1,N 
60      PRINT  19,    X(1,J),    >3A(1,J),    BA(2,J) 

PRINT  8,    SBA(l),    STA(2) 

DO  70   1=1,2 

GG(I)=sSV(I)*.r;WXX(I)-SWX(I)**2 

CD(I)=(SGV(I)*SWXX(I)-SGVX(I)*SWX(I))/CC(I) 

Gr'.(I)=(SV(I)*SGVX(I)-SG7(I)*SV7X(I))/GG(I) 

GF(I)sGD(I)f0.5*GS(I) 

AA(i)=(3U(I)*SXX(I)-SX(I)*SXU(I))/(SM(I)*SXX(I)-SX(I)*SX(I)) 

BB(I)=(SMCI)*SXU(I)-S:C(I)*SU(I))/(3M(I)*SXX(I)-3X(I)*SX(I)) 
70     AB(I)=AA(I)+0.5*B3(I) 

DTA1= ST9  CI ) + SRR  (  2  ) 

DTA2=3  ,T.(1  )+SVRO  ) 

DT32*SRR(1)+S7V(2) 

DTG1=S  V(1)*SUR(2) 

DT:2=SUR(1)+SUV(2) 

A12  (1  )=  (  DTG1  *D  TB2-DTC320TA2  )/(  DTA1  *DTB2-DTA2**2  ) 

A21  (1  )=  ( DTAl*0TC2-DTA2*DTCl )/(  <TA1*DT32-DTA2* *2  ) 

DO  80  J=1,N 
SQ=:S0^GALOZ(l,J)*  (2,J)**2*(l.-2.*X(l,J))+2.^5ALOE(2,J)*X(2,J)*(l,J)**2 
SW=SW+GAL0E(1,J)*2.  ^X(l,J)^:(2,J)--*2+GALOC(2,J)'X(l,J)*^2(l.-2#*X(2,J)) 

S0=80+X(2,J)**4*a.-2.>r(l,J))*^2+4.*X(l,J)*»4^X(2,J)**2 

SR=  ■'R*2.*(l.-2#*X(l,J))*xa,J)*X(2fJ)*(X(2tJ)   3-X(l,J)**3) 
80   SS=SS+4.*X(ltJ)**2*X(2,J)**4+X(l,J)**4*U.-2.*.:(2tJ))**2 

A12<2)«(SQ*SS~SR*SW)/(SO*SS-SR*SR) 

A2i(2)=(so*sw-s^*5:^)/(so*ss-r 

DO  90  1=1,2 

READ  16 

PRINT  30 

P".INT  16 

DO  85  J=1,N 

GAL0^(I,J)=(1.-X(I,J))**2*(AA(I)+BB(I)*X(I,J)) 

GBLOi-J(I,J)=(l.-X(I,J)**2*(3D(I)+CE(I)*X(I,J)) 

GCL0C(1,J)=(1.-X(1,J))**2*(A12(1)+2.*(A21(1)-A12<L))*X(1,J)) 

GC^00(2tJ)=(l.-X(2,J))**2*(A21(l)+2#*(A12(l)-A21(l))*X(2,J)) 

Gr>LOC(l,J)=(l.-X(l,J)y*2*(Al2(2)+2.*(A2](2)-Al2(2))*X(l,J)) 

G-)LOG(2,J)=a.-xr2,J))**2*(A21(2)+?.*(A12(2)-A21(2))*X(2,J)) 
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GEL0C(l,J)=(l.-X(l,J))**2*(SBA(l)+2.*(SBA(2)-S3A(l))*X(lfJ)) 
GELOC(2,JO=(1.-X(2,J))  *2*(SBA(2)+2.*(SBA(l)-SBA(2))*X(2fJ)) 
GAA(IfJ)=EXP(2.302585*GALOC(I,J)) 
GA3(I,J)=:EXP(2.3025n5*GBLOr;(I,J)) 
GAa(I,J)=EXP(2.302585*GCLOG(I,J)) 
GAD(I,J)=EXP(2.302535*GDLOC(I,J)) 
GA~(I,J)=EXP(2.302585*GELOC(I,J)) 

Sl(I)=Sl(I)+(C(GAA(I,J)-GA(I,J))/GA(I,J))**2)**0.5 
S2(I)=S2(I)+(((GAB(I,J)-GA(I,J))/GA(I,J))**2)**0.5 
S3U)=S3(I)+(((GAC(I,J)-GA(I,J))/GA(I,J))**2)**0.5 
S4(I)=S4(I)+(((GAD(I,J)-1A(I,J))/GA(I,J))**2)**0.5 
85  S5(I)=S5(I)  +  C((GAT^(I,J)-GA(ItJ))/GA(I,J>)**2)**0.5 
S1(I)=S1(I)/SM(I) 

S2(i)=S2(i)/.:.:(i) 
S3(I)=S3(I)/SM(I) 
S4(I)=S4(I)/SM(I) 
S5(I)=S5(I)/SM(I) 
PRINT  9 
11=1 
12=2 
13=3 
14=4 
15=5 

TEMP=AA*2) 
AA(2)=A3(2) 
AB(2)=TEMP 
tP=GD(2) 
CD(2)=GF(2) 
GF(2)=TEMP 

PRINT   7,    II,    AA(I),   AB(I)t    S1(I) 
PRIMT   7,    12,    31(1),    3?(I),    S2(l) 
PRINT  7,    13,   A12(1),   A21(l),    S3(I) 
BftWT   7,    14,    A12(2),    A21(2),    S4(I) 
90     PRITT   7,    15,    STlA(l),    SBA(2),    S5(I) 
GO  TO  10 
END 
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b) 

Title:  Program  for  the  calculation  of  the  ternary  constant  and 

ternary  activity  coefficients. 

Comput-r:   IBM  1620  with  50,000  positions  of  storage. 

Program  language:     Fortran. 

Description:  The  program  was  used  to  calculate  the  three  activity 

coefficients  in  a  ternary  mi::ture  by  Wohl*a  three-suffix  equations. 

Provision  has  been  ma  e  for  a  maximum  of  60  data  points.   The 
average  of  the  absolute  value  of  the  percent  deviation  between 
calculated  Tj  and  experimental  0   tov   each  component  was  calculated. 
If  the  ternary  constant  was  zero  in  the  input,  the  average  values 
of  the  ternary  constant  was  calculated  for  three  cases:  (1)  using 
all  data  points,  (2)  using  only  data  points  where  x  la  greater 
than  0.05,  and  (3)  using  only  data  points  where  xiX2 (1-2x3)  is 
greater  than  0.1,  and  x  is  greater  than  0.05.  The  calculation 
starts  after  setting  program  switch  1  in  the  en  position  and 
depressing  start  ke  . 

The  input  data  should  be  fed  into  machine  ir  the  following 
order: 

Card  1.   Comment  such  as  Benzene,  Gyclohexane,  Acetone  C=0 
less  than  49  characters  including  blanks  in  the  body  of  the  comment. 

Card  2.   A3.3,  A31,  A23,  A32,  A12,  A2i 

Card  3.   C,  qo/«H,  qsAl 

For  Margules  equation  set  o^/^i^q^/qi   =  1 

Having  program  switch  1  in  the  on  position  requires  that 
additional  input  data  be  fed  into  the  machine. 
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Card 

4.      MK,  K 

Card 

• 

5.     xu,  x2L,  yLi,  y2i 

* 

• 
• 
Card 

• 
• 

(4+K)   x1K,  x2K>  y1Kf  y2K 

MK  is 

!  a  fixed  point  constant  representing  the  form  of  exp- 

eriraental 

•  input.   K  is  also  a  fixed  point  constant  representing 

the  mimb^.i 

*  of  data  points 

After  all  x,  y  input  had  been  Pft«d|  the  following  cords  are 

fed  into  machine  depending  on  the  value  of  MK: 

for  MK=1 

Card 

(5+K)  Tl 

Card 

(6+K)   Alt  Blf  Ci 

Card 

(>+K)  A2,  B2,  C2 

Card 

(8+K)  A3,  33,  C3 

Card 

• 

(9+K)  Tx 

• 

• 
Card 

• 
(8+2K)  TK 

for  MK=2 

Card 

• 

• 

• 
• 

• 
(*+2K)  ilK,    ^K,  03K 

for  MK=3 

(5+K)  ln\Lf  ln&21,  lnX31 

• 

Card 

• 

• 
• 
Card 

• 

(4+2K)  lnOlK,  ln02K,  lnP3K 
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for  MK=4 

Card   (5+K)     log|u,   log»21,   log031 
•  • 

Card   (4+2K)  log\K,   logf2K,   IoS^3K 

All  output  is  presented  on  the  console  typewriter  in  the 
following  manner: 

Benzene  Cyclohexane  Acetone  C=0 

A13      A31      A23      A32       A12     A21 


Sum  of  A         Predicted  C 


Total         x  greater       u  greater 
than  0,05       than  0.1 
Data       Error   Data   Error     Data   Error 


Component  1 

Component  2 

Component  3 


•  •  •  •       ••••  •••• 

•  •  •  •      ••••  •••• 

•  •••      •  ♦  •  •  •••• 


If  C  is  zero  in  the  input  data  then  the  following  output 

is  typed, 

Total       x  greater       u  greater 
than  0.05       than  0.1 

Component  1 

component  *..   •**.   .»••   »•«•    . . . .     ...«    ••>. 

Component  3. 

Note:   Data  =  Number  of  data  points  / 

Error  =  Absolute  values  of  the  average  deviation  of  0   - 

from°exp« 
Sum  of  A  =  A]_3  +  A31  +  A23  ♦  A32  +  Ai2  +  A2]_ 

Predicted  C  •  0.5(A2^  -  A^2  +  A^3  -  X^   +  Ao->  -  A23> 
After  obtainirg  the  complete  output  for  a  system,  it  may  be 
desired  to  try  a  new  Margules  binary  constant  or  ternary  constant. 
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This  can  be  accomplished  through  setting  sense,  witch  L  in  the 
off  position,  x»rhich  allow  one  to  read  in  the  new  constant  without 
reading  x,  y,  and  ». 
Console  setting  for  operation: 

Overflow  check  switches  -  program 

All  other  check  switch  -  stop 

Sense  switch  1  -  on:  title  card  and  two  constant  cards  and 

x,  y,  ft  cards. 

off:  read  title  card  and  two  constant  cards 
only. 

Sense  switch  2,  3  and  *:■  -  immaterial. 


VI 


BLOCK     DIAGRAM      FOR     USE     WITH     IBM    1620 
FOR     COMPUTING     TERNARY     CONSTANTS    AND 
THE     ACTIVITY    COEFFICIENT 


START 


READ.TITLE 
COMMENT  ,A's 
C,  q2/qi,etc. 


CALCULATE      SUM      OF     A, 
PREDICTED     C,   B's 


PRINT       SUM     OF    A, 
PREDICTED     C 


SWV     OFF 


READ, MK  ,  K 
XI,X2,  Yl  ,Y2 


i 


CALCULATE    log  Vs 
BY    MARGULES    EQ. 


OFF 


•© 


6  0  0  0 
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READ, IT 

READ,Ai,Bt,Ci, 
A2£2,C2, 
A3>B3£3 

READ,  T 


CALCULATE 
EXPERIMENTAL 
"it  log  "Vt 


READ,7l,y2,>3 


CALCULATE 
EXPERIMENT, 
log"/,,  log  t2, 
log  i^__ 


READ,  InTfi, 
In  72,  In  K-j 


I 


CALCULATE     AVERAGE 
ABSOLUTE    DEVIATION 
BETWEEN    CALC.  1    AND 
EXP.   7    FOR     THREE 
CONCENTRATION     REGION 


PRINT    AVE. 
ABS.   VALUE 
OF    DEV. 


CALCULATE    TERNARY 

CONSTANT    IN    THREE 

DIFFERENT    CASES 

PRINT    AVE. 

ABS.  VALUE 

.OF   DEV.  AND 

TERNARY 

CONSTANT 

-1  START) 


A 


0  CALC^JLATION  OF  TERNARY  CONSTANT     ACTIVITY  COEFFICIENT 
DIMENSION  AA(3),AB(3),AC(3),T(60),P(3,60),GA(3,60),GALN(3,60) 
DIMENSION  A(3, 3), 0(3, 3), 3(5,3),      0)  ,Y(3,60)  ,GAl(3f  0)  ,  a(3,60) 
DIMENSION  GAL0E(3,  ))f  ,    0)  ,V(3,r>0)  ,SX1(3)  ,SK2(3)  ,SK3(3) 

DI      )N  S1(3),S2C3)/:3(3),U(3,60),TERNA(3,G0>,.%(3),S5(3),S6(3) 

ilON  G31(3,G0),GC1(3,60),G01(3,'0),G1(3,50) 
^1  .     I  SK4(3),  SK5(3),  SK6(3) 

1  FORMAT  (12,  12  ) 

2  FORMAT  (  EL2»6,  El  2.  6  ) 

3  FORMAT  (  E12.6,  12*6,  E12.6  ) 

h      FORMAT  (  E1.2.^,  E12.6,  E12.6,  E12.6  ) 

5  FORMAT  (9HC0MP0NENT  13,  E12.G,E12.6,^12. ' ,E12.6,E12.6,E12.6  ) 

6  FORMAT  (E12.6.R12.6,  EL2.6,  El 2 .  ,  12.6,  El?..   ) 

7  FORMAT  (/5X  3HA13  9X  3HA31  9X  3HA23  9X  3HA32  9X  3HA12  9X  3HA21) 

8  FORMAT  (49H  ) 

9  FORIiAT  (///) 

10  PAUSE 

IAD  8 
PRINT  8 
PRINT  7 

:AD  6,  A(1,3),A(3,1),A(2,3),A(2,2),A(1,2),A(2,1) 
PRINT  6,  A(1,3),A(3,1),A(2,3),A(3,2),A(1,2),A(2,1) 
RSA)  3,  CA,  ^(2,1),  o(3,l) 
IF  (  ^(2,1)  )  12,  11,  12 

11  0(2,1)=A(2,1)/A(1,2) 
0(3,1)=A(3,1)/A(1,3) 

12  0(1,1)=1. 
0(2,3)=o(2,l)/o(3,l) 
O(3,2)=o(3,l)/O(2,l) 
Q(l,2)=l./o(2,l) 
0<l,3)=l./o(3tl) 

SA2=A(2,1)+A(1,2)+A(3,2)+A(2,3)+A(1,3)+A(3,1) 
SAl=:(A(2,l)-A(V-)+-A(3t2)-A(2,3)+A(l,3)-A(3fl))*0.5 

13  FORMAT  (  12H   SUM  OF  A  12H  PREDICTED  C  ) 
PRINT  13 

PRINT  2,    SA2,    3A1 

15  FORMAT   (22 X  5HT0TAL   15a  12HX   BIGER   0.05   13X  11HU   BIGKft  0.1) 

16  FORMAT    (   16X4HDATA   8X5F^l^OR7X4HOATA8X5HERROR7X4HDATA8X5HERROR) 
DO  27   J=l,3 

L=J+1 

IF  9L-3)   22,    22,   21 

21  L=J+l-3 

22  M=J+2 

IF   (M-3)   24,    24,    23 

23  M=J+2-3 

2k      B(1,J)=:2.*(A(L,J)*0(J,L)-A(J,L)) 

B(2,J)=2.*(A(M,J)*Q(J,M)-A(J,M)) 

B(3,J)=A(JfM)+A(L,J)*0(J,L)-A(M,L)*0(J,M)-GA*0(J,l) 

B(4,J)=2.*(0A*1(J,1)+A(M,J)*^(J,M)-A(J,M)) 
27      3(5,J)=2.*(A(M,L)*0(J,M)-A(L,M)*n(J,D) 

IF    (    SENSE   SWITCH  1    )    28,    31 
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28  READ  1,  MK,  K 
00  29  N=0.,K 

m   4.  XKl,N)tXL(2#N)tYatH)tY<2,N) 

X1(3,N)=1.-  0.(1,N)-  .   ,  I) 

29  Y(3,N)=l.-Y(lfN)-Y(?,N) 

J  30  J=l,3 
00  30  N=1,K 

30  Zl(.7,N)=a(J,N)*0(.7,l)/(Xl(l,N)+Xl':2,N)*0(2,l)  +  'a(3tN)*Q(3,l)) 

31  00NTI  IDE 
DO  36  J=l,3 
">0  3*  N=1,K 
L=J+1 

IF  (L-3)  33,  33,  32 

32  L=L-3 

33  M=J+2 

IF  (  -3)  35,  35,  34 
3^  M=M-3 

35  GAl(J,N}=2aa,N)**2*(A(J,L)+Zl(J,N)*B(l,J)) 
OBl<JtN>aZlCM,N)**2*(A(J,M)+a(J>N)*a(2tJ)) 
Gal(J,N)=Zl(L,N)*Zl(M,N)*(B(3,J)+ZlCJ,N)*3(i*,J)+Zl(M,N)*B(5tJ)) 

aCJ,N)=GAl(J,N)+ciai(jjN)+GOi(jtN) 
Gl ( J, N)*EXP(2. 3025*5*001 < J, N)) 
U(J,N)=Zl(L,N)*Zl(MtN)*(l.-2.*Zl(J,n}) 

36  V(J,N)=(U(J,N)**2)**0.5 
DO  37  1=1,3 
SK1(I)=0. 

SK2(I)=0. 
SK3fI)=0. 
SK4(I)=  . 
SK5(I>=0, 
SK6(I)=0. 
S1(I)=0. 

:i>»o. 

S3(I)=0. 
S4(I)=->. 
S5(I)=\ 

37  S6(I)=0. 

IF  (       SWITCH  I  )  38,  49 

38  3     m 

GO  TO  (  39,  43,  45,  47  )  ,  MK 

39  READ  2,  PI 
HO  40  1=1,3 

40  READ  3,  AA(I),  AB(I),  AC(I) 
00  41  N=1,K 

41  READ  2,  TOO 
DO  42  1=1,3 
00  42  N=1,K 

P(I,N)=r.\'p(2.302585*(AA(I)-AB(I)/(AC(l,>+T(N)))) 
GA(I  ,N)sPi*Y(i  ,N)/P(I  ,N)/X1  (I , 

42  GAL0E(I,N)=L0G(GA(IfN)V2.302585 
GO  TO  49 

43  CONTINUE 
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DO  44  N=l,K 

READ  3,  GA(1,N),  GA(2,N),  GA(3,N) 

DO  44  1=1,3 

44  GALO^(I,N)=LOG(GA(I,N)')/2.302585 
GO  TO  49 

45  GONTIN' ' 
DO  hr.   N=l, 

READ  3,  GALN(1,N),  GALN(2,N),  GALN(3,N) 

DO  46  1=1,3 

GALOE(I ,N)=GALN(I ,N)/2.302535 

46  GA(I,N)=E      ,N(I,N)) 
GO  TO  49 

47  C 

DO  48  N=1,K 

REA)  3,  GAL0E(1,N),  G.\T,OE(2,N>,  GAL0^(3,N) 

DO  48  1=1,3 

48  GA(t,N)=  !  !  (2.302585*OALrR(l,flm 

49  CONTIN 

DO  53  1=1,3 
DO  53  N=1,K 

:I,N)=(((GA(I,N)-G1(I,N))/GA(I,N))**2)**0.5 
IF  (  Xl(I,N)-0.05  )  52,  52,  50 

50  S2(I)=S2(I)+ER^xO'i(I,N) 
SK2(I)=3K2(I)+L0 

IF  (  7(1,10-0.1  ^  52,  51,  51 

51  S3(l)«S3<I)+BPROR(t,N) 
SK3(i;=SK3(I)+1.0 

52  Sl(I)«Sl(I)*BRROa(I,H) 

53  SKL(I)=SK1(I)+1.0 
PRINT  9 

PRINT  15 
[HT  16 
D  54  1=1,3 
S1(I)=S1(I)/SK1(T) 
S2(T)=S2(I)/SK2(I) 
S3(I)«S3(I)/SK3(I) 

54  PRINT  5,  I,  SK1(I),  SKI),  SK2(I),  S2(I),  SK3(I),  S3U) 
IF  (CA)  10,  60,  10 

60  CONTINUE 

DO  64  1=1,3 
TO  CA   N=1,K 

TERNA(I,N)=(GALO-:(I,N)-GDl(l,N))/ua,N)*(-l.) 
IF  (    :a(l.N)-0.05  )  63,  63,  61 

61  S4(I)=S4(I)+-?^NA(I,N) 
SK4ri)=SK4(I)+1.0 

IF  (   V(I,N)-0.1  )  63,  62,  62 

62  S5(I)=S5(I)+TERNA(I,N) 
SK5(I)=SK5(I)+1.0 

63  S6(I)=S6(I)+TET>NA(I,N) 

64  SK6(I)=SK6(I)+1.0 
PRINT  9 

PRINT  65 


PRINT  15 

65  FORMAT  (  16HTTSRNARY  CONSTANT  ) 
00  66  1=1,3 
S4(I)=34(I)/SK4(I) 
S5(I)=S5(I)/SK5(T) 
S6(I)-S6<I)/SK6(I) 

66  print  5,  I,  SK6(I),S6(I),Sl^(I),Sft(I),SK5(I),S3(I) 

GO  TO  10 
END 
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ABSTRACT 

The  purpose  of  this  work  is  to  u     >hl*s  ternary  Margules 
three-suffix  equations  to  evaluate  the  activity  coefficients 
in  ternary  systems.   Binary  constants  for  the  equations  were 
obtained  from  the  three  constituent  binary  systems.   Using  the 
method  of  Least  squares,  several  methods  of  determining  binary 
constants  were  evaluated.   Inserting  these  six  independent  binary 
constants  into  the  Margules  equation  a  ternary  constant  waa 
evaluated  using  the  observed  ternary  activity  coefficients.  Eleven 
non-aqueous  ternary  systems,  found  in  the  literature,  were  used 
to  determine  the  need  for  a  ternary  constant. 

All  calculations  were  made  by  the  IBM  1620  computer  and  the 
computer  programs  for  determination  of  Margules  binary  constants 
and  calculation  of  ternary  activity  coefficients  are  presented. 

The  results  show  that  if  the  constituent  binary  systems  can 
be  well  represented  by  the  Margules  equation,  the  ternary  system, 
can  be  represented  using  the  three-suffix  Margules  equations 
without  a  ternary  constant,   -^he  precision  of  the  ore  -action  of 
activity  coefficients  in  the  ternary  system  is  about  the  same  as 
the  precision  with  which  binary  activity  coefficients  are  represent' 
ed  by  the  Margules  equation. 


